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Abstract 

When the number of particles N is finite, the noncohiding Brownian motion (BM) 
and the noncohiding squared Bessel process with index u > -1 (BESQ^'^)) are deter- 
minantal processes for arbitrary fixed initial configurations. In the present paper we 
prove that, if initial configurations are distributed with orthogonal symmetry, they are 
Pfaffian processes in the sense that any multitime correlation functions are expressed by 
Pfaffians. The 2x2 skew-symmetric matrix- valued correlation kernels of the Pfaffians 
processes are explicitly obtained by the equivalence between the noncolliding BM and 
an appropriate dilatation of a time reversal of the temporally inhomogeneous version of 
noncolliding BM with finite duration in which all particles start from the origin, N6o, 
and by the equivalence between the noncolliding BESQ^'^^ and that of the noncolliding 
squared generalized meander starting from N6o. 

Keywords Determinantal and Pfaffian processes, Eigenvalue distributions of ran- 
dom matrices, Noncolliding Brownian motion, Noncolliding squared Bessel process and 
generalized meander 

1 Introduction 

We consider one- dimensional particle systems called the noncolliding Brownian motion and 
the noncolliding squared Bessel process. The former is equivalent with Dyson's Brownian 
motion (BM) model with parameter P = 2, which was introduced as the eigenvalue process 
of the Hermitian-matrix-valued BM [S] corresponding to the Gaussian unitary ensemble 
(GUE) of random matrices [5^ [TT| . The latter is a one-parameter family indexed by z/ > 
—1 and abbreviated as noncolliding BESQ^*^-* [22]. When the number of particles is finite, 
N e N = {1,2,3,...} and u G No = NU {0}, the noncolliding BESQ^'') realizes the 
eigenvalue process of the matrix-valued diffusion process called the Laguerre process (or 
complex Wishart process) whose distribution at each time describes squares of singular 
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values of (A^ + u) x N random matrices in the chiral Gaussian unitary ensemble (cliGUE) 
[SnillS], and when z/ = 1/2 (resp. u = —1/2), it expresses the stochastic evolution [17j of 
the squares of positive eigenvalues of 2N x 2N random matrices in the Gaussian ensemble 
of class C (resp. class D) [H E]- (See [HlElESliaillHlllSlSSlEZlllliailOllZlE^ 
related interacting particle systems.) 

In the previous papers [20l [211 122], it was proved that if the number of particles is 
finite, G N, these two processes are determinantal for arbitrary initial configurations 
^(■) = J2f=i ^^ji') with Xi < X2 < ■ ■ ■ < xjy, Xj & A,l < j < N, where 6y{-) denotes the delta 
measure on y; 6y{x) = 6xy, A = M for the noncolliding BM and A = = {a; G M : x > 0} 
for the noncolliding BESQ'-'^^. Here, given a fixed initial configuration ^, a process is said to 
be determinantal, if there is a function K(s, x; t, y) such that it is continuous with respect 
to {x,y) for any fixed {s,t) G [0, oo)^ and any multitime correlation function is given by a 
determinant in the form 

pHtu 4!; • • • ; tM, ^) = det , [K(t^, ); t„, x^-\ (1.1) 

1 ^ J ^ i V 77], , 1 ^ ^ i V 7T, 

l<m,n<M 

M G N, < ti < ■ ■ ■ < tM < oo, 1 < A^^ < AT, 1 < m < M, where x^^^ = {x^^\ . . . , xj^^) 
denotes the points at which observation is performed at time tm, 1 < m < M [19j. The 
function IK is called the correlation kernel and it determines finite dimensional distributions 
of the process through (11.11) . For a configuration ^(■) = J2f=i^xj{-), shift by u; G C is 



denoted by r^.C(') = X]j=i ^Xj+w{-) and dilatation by factor c > is denoted by c o ^(■) = 
Ylj^=i ^cxj{-)- The correlation kernels for the noncolliding BM and the noncolliding BESQ^*^-* 
are respectively given by 

K^{s,x;t,y) = lim / du (p dzp{s,x\z)—^^^ p{—t,iu\y) 

elo 2m JM\[-s,e] JcuO tu- z 

-l{s>t)p{s-t,x\y), (x, y) GM^(s,t) G [0,oo)^ (1.2) 

K«(s,x;t,i/) = lim— [ du I dzp^''\s, x\z) ~ ^^ p('^)(-t, u\y) 

eio 2m /c4C) u- z 

-lis > t)p^''\s - t, x\y), (x, y) G (0, oof, [s, t) G [0, oof, v > -1, 

(1.3) 

where i = a/— 1, Cz'{C,) denotes a closed contour on the complex plane C encircling the 
points in supp ^ = {x G A : ^({x}) > 0} once in the positive direction but not the point z', p 
and p^^-* are the extended versions of transition probability densities of the one- dimensional 
standard BM [21] and the BESQ^'^) [22], 

I, J tGM\{0},x,yGC, 
V{t^y\x) = { ^2n\t\ (1.4) 

6{y-x), t = 0,x,yeC, 
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p^''\t,y\x) 



_L e-(^+^)/2*J. (^^^ , t G R \ {0}, X G C \ {0}, y G C, 

^ t G M\ {0},x = 0,1/ G C, (1.5) 



(2|t|)-+ir(z/ + i) 



, 5(y - x), t = 0,x,?/ G C, 

with the Gamma function T[z) = e~'^u^~^du and the modified Bessel function Iu{z) = 
Yl'^=oi^/'^y"^~^'^ /{^(^ + l)r(?2 + 1 + u)}, is an entire function having supp ^ as the zero 
set expressed by the following Weierstrass canonical product with genus 



n.w-n(i-i)- n (i-f)*". -c. (1.6) 

a'gg aiGsupp g 



and is the indicator of a condition oj; l{cj) = 1 if a; is satisfied and l(a;) = otherwise. 
In fll.Sp we have defined z" to be exp(z/log^), where the argument of z is given its principal 
value; = exp[z/{log |z| + A/^arg(z)}], — tt < arg(z) < n. We say that the correlation 
kernels, which are asymmetric, K(s, x; t, y) ^ ]K(t, y, s, x) for s 7^ t, as fll.2p and (11 .Sp . are 0/ 
Eynard-Mehta type ^M\U^^- 

For AT G N, e = Ef=i , a; = (xi, . . . , a;,v), /3 > 1, a > -1, a2 > 0, let 

-.-Af{/3(Af-l)+2}/2 

/^!v'l^(0 = e-l^l^/->^(a.)|^ (1.7) 

-N{l3{N-l)+2{a+l)} ^ 

= -i^^ l[{x^e~^^/"^')\h^{x)\^, (1.8) 

j=l 



where 

J) 

Kl.k<.l\ " 

i<j<k<N 

|ajp = ^27=1^]^ the normalization factors are given by 



/ijv(a;) = If (xk—Xj)= det [x^^" 

J-J- ^ ■'^ l<j,k<N ■> 



(^)_ (27r)^/^ -,^ r(j/3/2 + l) 
^ iV! IjL r(/3/2 + l) ' 

^^^^^ _ 2iv{/^(A^-i)+2(a+i)}/2 jv ^^^.^^2 + l)r(j/3/2 + a - /3/2 + 1) 



N\ J-jL r(/^/2 + 1) 



They are the probability density functions of random configurations S = ^Xj and 

S = Yl!j=i which the configuration spaces of particle positions X = (Xi, . . . , X^) and 

X = (Xi, . . . , Xtv) are given by 

X G = {a; = (xi, . . . , Xat) : Xi < X2 < ■ • • < Xn}, 

X G = {x = (xi, . . . , Xat) : < Xi < X2 < ■ ■ ■ < Xat}. 
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In particular, when /3 = 1, 2 and 4, fll.7p gives the distributions of eigenvalues of N x N 
Hermitian random matrices in the Gaussian orthogonal ensemble (GOE), GUE, and the 
Gaussian symplectic ensemble (GSE) with variances a^, respectively [321, [TT]. Similarly, for 
u e No, (|L8|) with (/3,a) = (1, (z/-l)/2), (2, z/), (4,2i/+l) give the distributions of squares of 
(distinct) singular values of (A^ + i>) x N random matrices in the chiral Gaussian orthogonal 
ensemble (chGOE), chGUE, and the chiral Gaussian symplectic ensemble (chGSE), respec- 
tively [Sni SHI HE]. Moreover [IB], (i) with a) = (l,z//2) was called 'the Laguerre 
ensemble (3 = 1 initial condition' in [12j, (ii) fll.8p with (/3,a) = (1,0) gives the distribution 
of squares of positive eigenvalues of 2N x 2N random matrices in the Gaussian ensem- 
ble of class CI studied by Altland and Zirnbauer [HE], (iii) (11.81) with (/3, a) = (1,-1/2) 
gives the distribution of squares of positive eigenvalues of 2N x 2N random matrices in the 
Gaussian ensemble of 'the real-component version of class D' of the Bogoliubov-de Gennes 
universality class [17j, and (iv) (11.81) with (/3,a) = (4,0) and (4,2) give the distributions of 
squares of distinct eigenvalues of 2N x 2N random matrices in the Gaussian ensembles of 
class Dlll-even [H [2] and of class Dlll-odd |13j|, respectively. We write the expectations of 
measurable functions of S and S, F{E) and F{E), with respect to distributions (II. 7p and 
(11.81) as E^|^2[-^(2)] and E^'"][F(S)], respectively. We can say that [T2] the distributions 
with (3 = 2 have unitary symmetry and those with (3 = 1 do orthogonal symmetry. 

Recently, we studied the noncoUiding BM and the noncoUiding BESQ^'"-', u > —1, starting 
not from any fixed configurations but from the distributions having unitary symmetry, /u)^ ^2 

and respectively. We showed that in these cases the determinantal structures of 

multitime correlation functions are maintained but the correlation kernels are replaced by 
the time shift t — )• t + o"^ of the correlation kernels for the special initial configuration 
^ = N60, i.e., the configuration in which all particles are put on the origin [14j. That is, 
the equalities 

p "•'^ {li, X^^, . . . , lM,X^^J — -Cjjy^^2 P {1^1, Xj^^, . . . , iM, ^Nm> 
l<m,n<M 



,(2,1-) 



det [Kf«(t„ + a^xf^^„ + a^x^'^))], z/ > -1, (1.9) 

l<m,n<M 

hold for any M G N, < ti < • • ■ < tM < 00, > 0, a;^^ G W^^ or a^J^^ G W+^, 
I < N.m < N,l < m < M. We should note that K^^" and K^^o are the correlation 
kernels known as the extended Hermite kernel and extended Laguerre kernel, respectively 

In the present paper, we report the cases when the noncoUiding BM and the noncoUiding 
BESQ*^*^^ start from the distributions having orthogonal symmetry, ^^^^2 and respec- 
tively. For G N and a skew-symmetric 2A^ x 2A^ matrix A = {ajk), the Pfaffian is defined 
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as 



1 >r-^' 

Pf (A) = Fii<j<k<2N{ajk) = J^Z^ sgn(7r)a^(i)^(2)a^(3)7r 



(4) ■ ■ ■ 0-Tr(2N-l)TT{2N), 



:i.io) 



where the summation is extended over all permutations vr of (1, 2, . . . , 2N) with restric- 
tion 7r(2/c — 1) < 7i{2k), k = 1,2, . . . , N. The main result of the present paper is the fact 
that for any cr^ > we can explicitly determine the 2x2 skew- symmetric matrix-valued 
correlation kernels 



A{s,X]t,y;a ) = 
{x, y) e (s, t) G [0, oof, and 

A'~'''^'>{s,x-t,y; a') = 



Au{s,x;t,y;a'^) Ai2(s, x; t, a^) 
-Ai2(t, y; s, x; ; a^) A22(s, x; t, y; ; a^) 



^[T\s, x; t, y; a^) A^g'"^ {s, x; t, y; a^) 
-AS2 (t, y\s,x; (T^) A^a'"^ a;; t, (x^) 



:i.l2) 



(x.y) G (0, 00)^, (s, t) G [0, 00)^, with k = 2{u — a) such that 



p^{ti,X 



(1). 



l<m,n<M 



'''N,a'^ 



■fnr -r^^^)^ 
) I'M; ■^^ATaj- > 



) I'M? •'^ATAf y 



:i.i3) 



Pfi<i<^^,i<fc<^n[A(''"^(tm,a;;."^^;t„,xl"^;a2)], > -l,a G (-l,z/], (1.14) 



l<m,n<M 



hold for any < ti < ■ ■ ■ < tM < 00, x^^^ G 



'Nr, 



or a:; 



(m) 



G W 



1 < A^™ ,< iV, 1 < 



m < M. As an analogue of a determinantal process, an interacting particle system is 
said to be a Pfajfian process, if any multitime correlation function is given by a Pfaffian 
[531iIlEnil3Iil37llI2lEHlE2[ISl[ISlE5]. Then we can state that noncolliding diffusion 
processes, which are determinantal processes if they start from fixed initial configurations, 
behave as Pfaffian processes when they start from distributions having orthogonal symmetry. 

The present paper is organized as follows. In Section 2 preliminaries and main results 
are given. In Section 3 is devoted to proofs of results. 



2 Preliminaries and Main Results 

Let N E N,0 < T < 00, and choose an initial configuration ^ = Sxj,Xi < X2 < ■ ■ ■ < 

xn- Then consider the A^-particle system of one- dimensional standard BMs starting from 
^ at time t = conditioned never to collide with each other during time period (0,T] [TB] . 
If the initial configuration is ^ = A^^o, that is, all A^ particles start from the origin, we can 
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show that the multitime joint probabihty density function for arbitrary M + 1 sequence of 
times < ti < ■ • • < < tM+i =T, M G N, is given by the formula 

M N 
X n f{tm+l-tm,x'^"'+''>\x'^"'y)hMix^''>)l[p{h,xf^\0), (2.1) 
m=l j=l 

where ^^""^ = ^J^^ = (xS"\ . . . , x57^) G W^, 1 < m < M + 1, 

fit,y\x)= det a;, G W^, t > 0, 

l<j,A;<A' 

and CN,T{t) = 7r^/2{nf=ir(j72)}-iT^(^-i)/4t-^(^-i)/2 The process, whose finite 

dimensional distributions are determined by the formula fl2.ip . is temporally inhomogeneous 
[T6| [23] . In this paper we call it 'the noncolliding BM with duration T starting from A^^o' 
and express it by (ST(t),t G [0,r],P^^o). 

When we take the limit T — > oo, we have a temporally homogeneous system [161 fT9] . 
which we simply call the noncolliding BM (starting from NSq). For the noncolliding BM, 
multitime joint probability density function is given by the following for an arbitrary initial 
configuration ^ with ^(M) = G N, 

M-l 

/(ti, . . . ; tM, e^^^) = hr.ix'-''^) J] /(Wi - ^n^, x^^'+'^\x'-^'>)h^^\h, x^'^- 0, (2.2) 



?Tl=l 



< ti < ■ ■ ■ < tM < C)0, with 



4+Ht,Z/;0 = ^ det 



1 / P(^,l/fck) 



27^^/c{f,) rixg^.l^-a;) 



, t>0,2/GW^, 



where, for a given initial configuration ^ = X^jLi^^jj^i — ^2 < ■■■ < a; at, we define 
~ Yli'k=i^^k^^ ^ j ^ N, and C(^) denotes a closed contour on the complex plane C 
encircling the points in supp ^ = {x G A : ^{{x}) > 0} once in the positive direction [5| 1^. 
The noncolliding BM starting from ^ is the temporally homogeneous process, whose finite 
dimensional distributions are determined by (12. 2p . and is denoted by (E(t),t G [0, 00), P^) in 
this paper. We can prove that H(t, ■) = X^^i (<)(■)> ^ — solves the following system of 
stochastic differential equations (SDEs), 

dX,{t) = dB,{t) + x (t)-X (tV 1^^'^^' ^^0' 

with independent one-dimensional standard BMs {Bj(t)}^^^, which is the (3 = 2 case of 
Dyson's BM model [T9l[23]. 
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In [18], a temporally inhomogeneous noncoUiding diffusion process was introduced, which 
is called the noncoUiding squared generalized meander with duration T. It is a two-parameter 
family of processes indexed by z/ > — 1 and k G [0, 2(z/ + 1)) starting from the configuration 
NSq, N ^ N, which includes the processes studied in [12] and [34] as special cases. This family 
of processes is denoted here by {'E^''^\t),t G [0, T], P^''°). The multitime joint probability 
density function is given by 



• • ■■^tM,C'^''^;tM^i,e'^''^) = <^ (ti)sgn(/.^(:.(^'^+i))) 



M 



TV 



(2.3) 



m=l 



< ti < ■ ■ ■ < tM < tM+1 = T, where 

f 1 /7/\('^-«)/2 



yV-K/2^-y/2t 

, 5(2/ -a;). 



t > 0,x > 0,?/ > 0, 

t > 0,a; = 0,2/ > 0, (2.4) 

t = 0, a;, ?/ > 0, 



f^'^KtM^) = det [p(^''')(t,2/,|a;fc)], a?,?/ e W+,t > 0, 

l<j,A;<TV 



and 



^ ^ 2n(Ar+K-l)Ar/2^-(TV-l)TV ^ 



r(i. + i)r(i/2) 



Lr(j72)r((j + l + 2z/-/s:)/2)- 



On the other hand, for the noncoUiding BESQ*^'^^ which is obtained as temporally homoge- 
neous limit T — 7- oo of S^^''^\t) [TTJ [18], the multitime joint probability density function is 
obtained for an arbitrary initial configuration ^ = Ylj=i ^Xj,^ ^ Xi < X2 < ■ ■ ■ < xn, N E N 
as 



M-l 



Pi{tu^^'\ . . . ; tM, e^^')) = /^Tv(a;(*-^)) n /^'^HWi - t^, a.("+^)|a;(-))/.!;'+)(ti, x^'^; 0, (2.5) 



m=l 



<ti < ■ ■■ <tM < OO, with [H [22] 



/('^)(t,?/|a;)= det [p(^)(t,y,|a;fc)], a;, G W+ , t > 0, 

l<j,fc<TV 



/^!;''"Ht,2/;0 = ^ det 

l<j,fc<TV 



1 

27Ti 



dz 



We write the noncoUiding BESQ'^'^^ starting from ^ as {'E!-''\t),t G [0,oo),P^). If we set 
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■) = ^f^^ X{t) = (Xi(t), . . . satisfies the SDEs 



dXj{t) = 2\J Xj{t)dBj{t) + 2{v + l)dt 

+AXj{t) ^ ^ , 1 < J < iV, t > 0, 

where {Bj{t)}jLi are independent one-dimensional standard BMs and, if — 1 < z/ < 0, the 
refiection boundary condition is assumed at the origin [22]. 

We write the noncolhding BM starting from the GOE eigenvalue distribution yuj^^^2 as 
_ (1) 

G [0, oo),P '^•'^'^) and the noncolliding BESQ'-'^^, u > —1 from the distribution having 

orthogonal symmetry, /i]y'"2,a G (— l,z/], as {E^'^\t),t E [0, oo), P'^^-'^^). 

In general, two processes having the same state space are said to be equivalent if they 
have the same finite-dimensional distributions, that is, if, for any finite sequence of times 
< ti < ■ ■ ■ < tM < oo, M E N, the multitime joint probability density functions coincide 
with each other [Hj. The key lemma of the present study is the following equivalence. 

Lemma 2.1 For > 0, let 

2 

c^2(t) = ^— , tG[0,oo). (2.6) 

Then 

G [0,oo),P<'^^) = (^-l-^oS.2(aV(t)),t e [0,oo),P^^»^ . (2.7) 

// the relation 

a = u - z/>-l, aG(-l,z/], (2.8) 

is satisfied, then 

(S('^)(t),t G [0,oo),P^-.'^^) = (^—L^ oS(r^(aV^(t)),t G [0,oo),P^^°^ . (2.9) 



Remark that 

0<s<t<oo >(r^c„2{s) > a^c^2{t) > 0. (2.10) 

Therefore, the RHS of (12. 7p and (12.91) are time reverses of the processes St(-) and 

with duration T = a^, followed by dilatation with factors l/Co.2(-) and l/Co-2(-)^, respectively. 

In I3S1II5] and [IB], it was proved that (^^(t), t G [0, T], P^^") and {Ep'"\t),t G [0, T], P^-^"), 
< T < oo, u > —1, K G [0,2(i/-|- 1)) are Pfaffian processes, respectively. Then by the 
equivalence (12. 7p and (12. 9 p of Lemma 12. ![ the following main Theorems are obtained. (For 
simplicity of expressions, we show the elements of the matrix- valued correlation kernels only 
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for the case that the number of particles is even. See, for example, |35] for the general 
theory of Pfaffian expressions of correlation functions.) Let if„ and LJ^ be the Hermite 
polynomial of degree n and the Laguerre polynomials of degree n with index v\ 



V21 



HJx 



n 



{2x 



k=0 



T{k + iy + l){n-k)\kV 



where [r] denotes the greatest integer not greater than r. For n G Z = {...,— 1,0,1, 2,...} 
and a G M we define 



n + a 
n 



' r(n + a + l) 

r(n + l)r(a + l)^ 
(-l)"r(-a) 



if n G N, a ^ Z_, 
if n G N, n + « G Z_, 



r(n + l)r(-n-a) 

0, if n G N, a G Z_, n + a G No, 

1, ifn = 0, 
0, if n G Z_, 



where Z_ = Z\No. 



Theorem 2.2 The noncolliding BM with a finite number of particles G N, starting from 

(1) 

the GOE eigenvalue distribution with variance cr^ > 0, (S(t),t G [0, oo), P'^^-'^^), is a Pfaffian 
process. When N is even, the elements of the matrix-valued correlation kernel U.ll\) are 
given by 



N/2-1 



k=0 



Au{s,x;t,y;(7'^) = ^ , B2k{s,x;a'^)B2k+i{t,y;a'^) - B2k+i{s,x;a'^)B2k{t,y;a'^ 



Ai2{s,x;t,y;cr'^) 

N/2-1 ^ 



k=0 



E 

k=N/2 



B2k+iis,x; a^)C2k{t,y; a^) - B2k{s,x; a'^)C2k+iit,y; 



B2k+i{s, x; cr'^)C2k{t, y; a^) - B2k{s, x; a^)C2k+i{t, y; 



ifs<t, 
if s>t, 



A22{s,x]t,y]a'^) = , , C2k{s,x;a^)C2k+i{t,y]a'^) - C2k+i{s,x;a'^)C2k{t,y;a'^ 

^ — ^ ai.(cr ' 



k=N/2 



(2.11) 
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with 



B. 



2k+i{s,x; CT^) 



2a^T{k + l/2)T{k + l), 



a^ + 2s 



2k 



X 



4a2 ) 

X { H2k+1 



X 



-x2/2{o-2+s) 



H2k-1 



Vcr^ + 2sJ a^ + 2s Wcj^ + 2s 



X 



X 



A;! 



{2i+l)\ \a^ + 2s 



a^ + 2s 



(2f+l)/2 



X 



vV2 + 2s 



_2-2fc+i 



0" 



a" 



a^ + 2s \a^ + 2s 



X 



gxV2(<T2+.)-xV{a2+2.)^_ 



2fc 



+ 2s 



k e Nn. 



(2.12) 



Theorem 2.3 T/ie noncolliding BESQ^'^\u > —1 wzt/i a finite number of particles N E 
N, starting from the distribution having orthogonal symmetry, /ij^'"2,cr^ > 0,a G (— l,z/], 

(l,a) 

{'E^'^\t),t G [0, oo), P^^.-^^), zs a Pfaffian process. Let k = 2{v — a). Then, when N is even, 
the elements of the matrix-valued correlation kernel M.12\] are given by 



Atr\s,x-t,y-a') 

N/2-l 



A^ii''\s,X]t,y]G'^) 

N/2-1 



Bir\s, X- a^)B^Zi{t, y; a^) - B^^ii^, x; a^)B^T'it, ^ 



2\ t}{v,k) 



2\ T}{^,n), 



fc=0 



S 4(U ^2lfi(^,^;^')c^r^(^,2/;^^)-5^rH^,^;cr^)c^£v1(^.2/;^ 



> t)p^^'''\s,x]t,y]a'^), 
A\^^'^\s,x]t,y]a^) 



k=N/2 "fc \" ) 



Cir\s, x; y; a^) - C'2ul[{s. x; a')C',r\t, y; a^) 



^2\M'^,i^) 



(u,k) 



2\^(i^,k)/ 



(2.13) 



10 



with 



-2. _2«(2A;)!r(2A; + 2 + 2z/-«;) 



A;! 



r(i/ + i)2 
1 



2^+^V{v + 1) (72(«+l)((72 + sy-" 



X 



J=0 



+ 2s 



(2A;)!r(2t/ - K + 1) 2 (f^^ + s)''"" /2A; + 2i/ - k + 1 
2^-ir(z/+ 1) {a^ + 2sY+\ 2fc + l 



^ gX/2(a2+s)-x/(o-2+2s)^, 



^2fe+i(*'^;'^^) — 



{2k + l)!r(2z/ - K + 1) 2 (t^^ + sY''" + 2V-K + 1 
2^-ir(z/ + 1) ^ ((t2 + 25)^^+1 V 2A; + 1 



X e 



a;/2((T2+s)_a;/(o-2+2s) 



j=2k 



2k--, 



r(j + 1) 



T{j + l + iy) ya'^ + 2s) ^a"^ + 2s J 



X \ 



(2.14) 



k e No, where 



^j,2k 



k — j + U — K 

k-j 

k + 2i/ — n fk — 2 — j + u — k\ /k — j + u — k, 
k-2-j )~\ k-j 



if k is even, 
, if k is odd, 



^]-2k-u + K-2 

[0-+l)/2] 



j > '2k, 



j,2k+i^- E bC^k + 3,2£-l) 



with 



e=k+i 

( (n-ra)l2 



j-2i-U + K-l 

j-2e + i 



, j>2k + l. 



h{m, n) — < 



m + 2i + 2u-K A ^ 

II , ij m, m are odd and m < n. 



m + 2i 



1, 
I 0, 



if m, n are odd and m > n, 
otherwise. 



and 



{v,k) / J. 2\ 



2 I \ U—K 

+ s 



a^ + t 
a^ + t 



^,/2(.2+.)-,/2(.2+t)^(.,K)(^ _ ^^ ^1^)^ ^ ^ 



-(•^-1^) / „2 



0" + S 
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for s > t,y > 0. 



3 Proofs of Theorems 
3.1 Proof of Lemma 12.11 

When the initial configuration ^ = ^f^i^x^ is distributed according to yuj^^^2,o"^ > 0, there 
is no multiple point in x = {xi, . . . , x^) with probability one, i.e., P^^'^ai^ G W^] = 1. In 
this case h^j^\t,y; C,) = f{t,y\x)/hN{x), and we can confirm the equality 

N 

x/(o-^ - a^c<^2(ti),a;|c^2(ti)a;(^))sgn(/iAr(a3)), 

where, for c > 0, t/ = . . . , un) G W;^, we put cy = {cyi, . . . , ci/n)- Similarly, we can see 
the equalities for 1 < m < M — 1, 

f (tm+l tfYi) x^ ^ I aj^ ^ ) 

x/(a2c,2(tj - a2c,2(t„+i),c,2(t^)a;('")|c,2(t^+i)a;(™+i)). 

Then for any M G N, < ti < < ■ ■ • < < oo, ^ = T.f=i^x,,x G W^, and 
^(m) ^ ^^^^^ a;(™) G W^, 1 < m < M, the equality 

/i!i!.2(0/(ti,e(^);...;tM,e(*^)) 

M 
m=l 

(3.1) 

holds. Integration of the LHS of (13.11) over x G W;^ gives 

/ cia./.Sl2(0/(^i,e^'^;---;Uf,^(^')) = EW^Jp^(ti,e(^);...;tM,e(*'))" 

= /-^(ti,e(^);...;tM,e(^)), 

and that of the RHS of (13. ip gives the multitime joint probability density function of the non- 
colliding BM starting from A^^q with duration T = cr^, in which observations are performed 
at M times in the reversed order < a'^c„2{tM) < o'^c^'^itM-i) < ■ ■ ■ < a^c^2(ti) < cr^, 
multiplied by the scale factors 0^2 (tm), 1 < m < M. Then the equivalence of the processes 
(12.71) is concluded. In a similar way, we can prove (12.91) . 1 
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3.2 Proof of Theorem [2721 

For a sequence (A^m)m=i of positive integers less than or equal to N, the {Ni, . . . , Nm)- 
multitime correlation function at M times < ti < • • • < < T of (Sr(t),t G [0,T],F^^°) 
is obtained from (12. ip by 



fjrp yii^ "^Ni 1 • • • 1 ''M? 



a; 



AT 



M • 



In [Mills], the functions Aii(s, x; t, ?/; T, ti), Ai2(s, a;; t, T, ti), A22(s, x; t, ?/; T, ti), < s, t < 
T, (x, y) G are given such that 

^N5o(. (1). 



Pfl<j<Ar^,l<fc<Ar„ 
l<m,n<Af 



By direct calculation we have found that 

Ajfc(o-^c^2 (s), c<^2 cr^c^2 (t), c^2 c^2 (tM)) 

= c,2 (s)-^/2c,2 (t)-i/2A^.fe(s, x; t, a^), (3.2) 

for (j, fc) = (1, 1), (1, 2), (2, 2), where A.^, (j, A;) = (1, 1), (1, 2), (2, 2), are given by i^JB) with 
(12112). Note that by definition of Pfafiian ffTTU]) . with any set of factors w^-, 1 < j < 2iV, 

2N 

Pil<j<k<2N{VjajkVk) = Y\_^j ^ Pf l<j<A:<2Af(ajfc)- (3.3) 

Then by the equality (12. 7p of Lemma HH] (see also (13.10 ). 



-^ATii • • • 1 I'M? •^iVJ^,^ , 



Pf 



l<j<Nm,l<k<Nn 

l<m,n<M 



A/ 



X 

^=1 



= Pfl<i<7V^,l<fc<7V„[A(tm, 4"''; ^fc"^! ^^)]' 
l<m,n<M 

and the proof is completed, i 



13 



3.3 Proof of Theorem [2731 

In [IB], the functions PS^{'^\s,X]t,y]T,ti), A^^2'^\s,x;t,y;T,ti), A'^2'^\s,x;t,y;T,ti), < 
s,t <T, (x, y) G (0, oo)^ are given such that 



^11 y^m^-'^j -i^n-i-^i^ i^^^l) -'^12 y^mi-'^j i^ni-^]^ i ? ■'1 
42 



,(1). 



'Ph<j<Nm,l<k<N.a 

l<m,n<M 



JT-iT K'^ni-^k I'^mi-'^j -'^22 V-m; •'^j ■j'^ni-^k ) 5 ''ly 



By direct calculation we have found that 

= c^2 (s)"^c<^2(t)~^A5^'''^(s, x; t, y; a^), 



(3.4) 



for {j,k) = (1,1), (1,2), (2, 2), where A^/^ (j, fc) = (1, 1), (1, 2), (2, 2), are given by 

with (12.141) . Then by the equality (12. 9 p of Lemma [2m and the property of Pfaffian (13.31) . the 

theorem is proved, i 



Acknowledgements The present author would like to thank T. Imamura for useful comments 
on the manuscript. A part of the present work was done during the participation of the present 
author in the ESI program "Combinatorics and Statistical Physics" (March and May in 2008). The 
author expresses his gratitude for hospitality of the Erwin Schrodinger Institute (ESI) in Vienna and 
for well-organization of the program by M. Drmota and C. Krattenthaler. This work is supported 
in part by the Grant-in-Aid for Scientific Research (C) (No. 21540397) of Japan Society for the 
Promotion of Science. 



References 

[1] Altland, A., Zirnbauer, M. R.: Random matrix theory of a chaotic Andreev quantum dot. 
Phys. Rev. Lett. 76, 3420-3424 (1996) 

[2] Altland, A., Zirnbauer, M. R.: Nonstandard symmetry classes in mesoscopic normal- 
superconducting hybrid structure. Phys. Rev. B 55, 1142-1161 (1997) 

[3] Bleher, P. M., Kuijlaars, A. B.: Integral representations for multiple Hermite and multiple 
Laguerre polynomials. Ann. Inst. Fourier. 55, 2001-2014 (2005) 

[4] Borodin, A., Ferrari, P. L., Prahfer, M., Sasamoto, T., Warren, J.: Maximum of Dyson 
Brownian motion and non-colliding systems with a boundary. Elect. Comm. Probab. 14, 486- 
494 (2009) 

[5] Borodin, A., Rains, E. M.: Eynard-Mehta theorem, Schur process, and their pfaffian analogs. 
J. Stat. Phys., 121, 291-317 (2005) 

[6] Brezin, E., Hikami, S., Larkin, A. I.: Level statistics inside the vortex of a superconductor and 
symplectic random-matrix theory in an external source. Phys. Rev. B 60, 3589-3602 (1999) 



14 



Delvaux, S., Kuijlaars, A. B., Zhang, L.: Critical behavior of non-intersecting Brownian mo- 
tions at a tacnode. Commun. Pure Appl. Math. 64, 1305-1383 (2011) 

Desrosiers, P., Forrester, P. J.: A note on biorthogonal ensembles. J. Approx. Theory 152, 
167-187 (2008) 

Dyson, F. J.: A Brownian-motion model for the eigenvalues of a random matrix. J. Math. 
Phys. 3, 1191-1198 (1962) 

Eynard, B., Mehta, M. L.: Matrices coupled in a chain: I. Eigenvalue correlations. J. Phys. A 
31, 4449-4456 (1998) 

Forrester, P. J.: Log-gases and Random Matrices. London Mathematical Society Monographs, 
Princeton University Press, Princeton (2010) 

Forrester, P.J., Nagao, T., Honner, G.: Correlations for the orthogonal-unitary and symplectic- 
unitary transitions at the hard and soft edges. Nucl. Phys. B 553 [PM], 601-643 (1999) 



Ivanov, D. A.: Random-matrix ensembles in p-wave vortices. (2001); arXiv:cond-mat/0103089 



Katori, M.: Characteristic polynomials of random matrices and noncolliding diffusion pro- 
cesses. arXiv:math. PR/1102.4655 

Katori, M., Nagao, T., Tanemura, H.: Infinite systems of non-colliding Brownian particles. 
Adv. Stud, in Pure Math. 39, 283-306, "Stochastic Analysis on Large Scale Interacting Sys- 
tems", Mathematical Society of Japan, Tokyo (2004) 

Katori, M., Tanemura, H.: Scaling limit of vicious walks and two-matrix model, Phys. Rev. E 
66, 011105/1-12 (2002) 

Katori, M., Tanemura, H.: Symmetry of matrix- valued stochastic processes and noncolliding 
diffusion particle systems. J. Math. Phys. 45, 3058-3085 (2004) 

Katori, M., Tanemura, H.: Infinite systems of noncolliding generalized meanders and Riemann- 
Liouville differintegrals. Probab. Theory Relat. Fields 138, 113-156 (2007) 

Katori, M., Tanemura, H.: Noncolliding Brownian motion and determinantal processes. J. 
Stat. Phys. 129, 1233-1277 (2007) 

Katori, M., Tanemura, H.: Zeros of Airy function and relaxation process. J. Stat. Phys. 136, 
1177-1204 (2009) 

Katori, M., Tanemura, H.: Non-equilibrium dynamics of Dyson's model with an infinite num- 
ber of particles. Commun. Math. Phys. 293, 469-497 (2010) 

Katori, M., Tanemura, H.: Noncolliding squared Bessel processes. J. Stat. Phys. 142, 592-615 
(2011) 

Katori, M., Tanemura, H.: Noncolliding processes, matrix-valued processes and determinantal 
processes, to be published in Sugaku Expositions (AMS); arXiv:math. PR/1005. 0533 



15 



[24] Katori, M., Tanemura, H.: Complex Brownian motion representation of the Dyson model. 
(2010); arXiv:math. PR/1008.2821 

[25] Katori, M., Tanemura, H., Nagao, T., Komatsuda, N.: Vicious walk with a wall, noncoUiding 
meanders, chiral and Bogoliubov-de Gennes random matrices. Phys. Rev. E 68, 021112/1-16 
(2003) 

[26] Konig, W., O'Connell, N.: Eigenvalues of the Laguerre process as non-colliding squared Bessel 
process. Elec. Comm. Probab. 6, 107-114 (2001) 

[27] Kuijlaars, A. B., Martmez-Finkelshtein, A., Wielonsky, F.: Non-intersecting squared Bessel 
paths and multiple orthogonal polynomials for modified Bessel weight. Commun. Math. Phys. 
286, 217-275 (2009) 

[28] Kuijlaars, A. B., Martmez-Finkelshtein, A., Wielonsky, F.: Non-intersecting squared Bessel 
paths: critical time and double scaling limit. arXiv:math.CA/1011.1278 

[29] Levin, B. Ya.: Lectures on Entire Functions. Translations of Mathematical Monographs, 150, 
Providence R. I.: Amer. Math. Soc. (1996) 

[30] Mahoux, G., Mehta, M. L.: A method of integration over matrix variables IV. J. Physique I 
(France) 1, 1093-1108 (1991) 

[31] Mahoux, G., Mchta, M. L., Normand, J. M.: Matrices coupled in a chain: II. Spacing functions. 
J. Phys. A 31, 4457-4464 (1998) 

[32] Mehta, M. L.: Random Matrices. 3rd edn. Elsevier, Amsterdam (2004) 

[33] Mehta, M. L., A. Pandey, A.: On some Gaussian ensemble of Hermitian matrices. J. Phys. A: 
Math. Gen. 16, 2655-2684 (1983) 

[34] Nagao, T.: Dynamical correlations for vicious random walk with a wall. Nucl. Phys. B658 
[FS], 373-396 (2003) 

[35] Nagao, T.: Pfaffian expressions for random matrix correlation functions. J. Stat. Phys. 129, 
1137 (2007) 

[36] Nagao, T., Forrester, P. J.: Multilevel dynamical correlation function for Dyson's Brownian 
motion model of random matrices. Phys. Lett. A247, 42-46 (1998) 

[37] Nagao, T., Forrester, P. J.: Quaternion determinant expressions for multilevel dynamical 
correlation functions of parametric random matrices. Nucl. Phys. B563[PM], 547-572 (1999) 

[38] Nagao, T., Katori, M., Tanemura, H.: Dynamical correlations among vicious random walkers. 
Phys. Lett. A 307, 29-35 (2003) 

[39] Noguchi, J.: Introduction to Complex Analysis. Translations of Mathematical Monographs, 
168, Providence R. I., Amer. Math. Soc. (1998) 

[40] Olshanski, G.: Laguerre and Meixner symmetric functions, and infinite-dimensional diffusion 
processes. J. Math. Sci. (New York) 174, 41-57 (2011) 



16 



[41] Pandey, A., Mehta, M. L.: Gaussian ensembles of random Hermitian intermediate between 
orthogonal and unitary ones. Commun. Math. Phys. 87, 449-468 (1983) 

[42] Rambeau, J., Schehr, G.: Extremal statistics of curved growing interfaces in 1+1 dimensions. 
Europhys. Lett. 91, 60006/1-6 (2010) 

[43] Rambeau, J., Schehr, G: Distribution of the time at which N vicious walkers reach their 
maximal height. Phys. Rev. E 83, 061146/1-27 (2011) 

[44] Revuz, D., Yor, M.: Continuous Martingales and Brownian Motion. 3rd edn. Springer, Now 
York (1998) 

[45] Schehr, G., Majumdar, S. N., Comtet, A., Randon-Furling, J.: Exact distribution of the 
maximal height of p vicious walkers. Phys. Rev. Lett. 101, 150601/1-4 (2008) 

[46] Sener, M. K., Verbaarschot, J.J. M.: Universality in chiral random matrix theory at ^ = 1 
and /3 = 4. Phys. Rev. Lett. 81, 248-251 (1998) 

[47] Tracy, C. A., Widom, H.: Differential equations for Dyson processes. Commun. Math. Phys. 
252, 7-41 (2004) 

[48] Tracy, C. A., Widom, H.: Nonintersecting Brownian excursions. Ann. Appl. Probab. 17, 953- 
979 (2007) 

[49] Verbaarschot, J.: The spectrum of the Dirac operator near zero virtuality for Nc = 2 and 
chiral random matrix theory. Nucl. Phys. B 426 [FS], 559-574 (1994) 

[50] Verbaarschot, J. J. M., Zahed, L: Spectral density of the QCD Dirac operator near zero 
virtuality. Phys. Rev. Lett. 70, 3852-3855 (1993) 



17 



